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ABSTRACT 

We review the pNGB quintessence models, and point out that the reason why the large decay 
constants f a > 0{l)M P i are really needed is to tame a tachyonic instability present for a 
wide range of initial vevs. Starting very close to potential maxima does not help because 
quantum fluctuations during early inflation at a scale Hi perturb the quintessence vev, 
displacing it from the maxima. This issue is quite interesting for pNGB dark energy in light 
of the recently discussed difficulties with embedding models with f a > Mpi in fundamental 
theory. A possible way around is provided by models with several ultralight pNGBs, which 
can drive a short burst of very late inflation together even if all of their decay constants 
obey f a < Mp\. Starting with their vev ~ f a , the pNGBs will hold each other up on the 
potential for a longer time period. Their effective dynamics is captured by a collective mode, 
containing admixtures of all of the rolling pNGBs, which behaves as an ultralight field with 
f e ff > Mpi. We point out that there may be potentially observable large scale disturbances 
in the sea of dark energy in such models. 
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"The time has come, " the Walrus said, 
"To talk of many things: 
Of shoes-and ships-and sealing-wax- 
Of cabbages-and kings- 
And why the sea is boiling hot- 
And whether pigs have wings. " 

(Lewis Carroll) 

There is mounting evidence that our universe may be accelerating [1, 2]. If gravity obeys 
Einstein's General Relativity (GR), this means that the universe is currently dominated 
by a dark energy component, contributing about ~ 70% to the critical energy density. 
To fuel cosmic acceleration dark energy ought to have negative pressure, satisfying roughly 
w = pi ' p £j —2/3 to fit the observations [3, 4], if it is the only agent responsible for supernovae 
dimming 1 . Because dark energy would control both the present state of the universe and the 
course of its future evolution, it is extremely interesting and important to try to elucidate 
the microscopic nature of dark energy. 

The first step in this quest is to devise benchmark models which are useful to compare 
with observations. Dark energy is usually modelled as a cosmological constant or as a time- 
dependent quintessence field [7, 8, 9, 10]. Each must be finely tuned to fit the observations 
[11]. Even if we leave aside the mystery of tuning the vacuum energy to nearly zero, and 
assume that the current acceleration is due to a quintessence field, it is hard to build realistic 
models. Quintessence requires more fine tunings, over and above the scale of dark energy 
density: a tiny mass itlq smaller than the current Hubble parameter Hq ~ 10~ 33 eV as well 
as sub-gravitational couplings to normal matter, to satisfy the Solar system gravity bounds 
[12]. This leaves but a few models which appear natural from the point of view of 4D 
effective field theory, in the sense that their masses and their couplings to the visible matter 
remain small in spite of the corrections arising in the loop expansion of quantum field theory. 
In these models quintessence is a pseudo-Nambu-Goldstone boson (pNGB) [13, 14, 15, 16], 
which in perturbation theory has only naturally weak, derivative couplings to visible matter, 
protected by a shift symmetry — > + C. Its potential, V = ji A [1 — cos(0// a )] + . . ., and 
specifically the mass term uiq ~ /i 2 / f a are generated by non-perturbative effects breaking 
the shift symmetry to a discrete subgroup, and hence are radiatively stable [13, 14, 15, 16]. 

It is often quoted that to serve as quintessence, the pNGB models need a large decay 
constant f a > M Pt [13, 14]. The idea is that the dynamics of quintessence should be akin 
to chaotic inflation [17, 18, 19, 20] where a dynamical field impersonates dark energy if the 
potential admits a slow roll regime. However, for models where f a > 0(l)Mp t there will 
generically be higher order corrections to the pNGB potential, coming as higher harmonics 
in the expansion over the instanton contributions, which typically spoil the flatness of the 
potential [20, 21]. Although one may invent models where such large values of f a appear 
and the potential remains flat [20, 22, 23], it has been questioned whether they can come 
from fundamental theory, e.g string theory [21]. Even though higher order corrections might 
still allow an e-fold or few of acceleration to unravel, which could be enough to account for 
observations, this poses some challenges for model building. 

1 See [5, 6] for a recent analysis of the observational status of some alternatives. 
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What of pNGBs with smaller f a 7 One might be tempted to try to squeeze some ac- 
celeration from them by arranging for the pNGB to dwell on, or near, a maximum of its 
non-perturbatively generated potential. Of course this entails very special initial conditions 
for the pNGB vev, namely its initial location on the potential peak. Also, one may try to use 
additional, explicit PQ symmetry breaking terms to shift the vacuum energy, changing the 
Lagrangian as in [24], which starts with the QCD axion and chooses f a ~ 10 12 GeV to argue 
that the axion may be dark matter. However, in this case dark energy is really an explicit 
cosmological constant term, as opposed to a dynamical quintessence, and we have nothing 
new to say about it. Note also that a pNGB with a mass m a ~ 10 _16 eV might weakly mix 
with photons, inducing a conversion of light from distant supernovae in the extragalactic 
magnetic fields, and thus also affect the observed supernovae dimming [25]. 

In this article we will review the dynamical constraints on the decay width f a needed 
for pNGB to behave as quintessence. We will demonstrate that the simple estimates of f a 
based on slow roll approximations are incomplete and somewhat misleading. The point is 
that the slow roll can set in only when the field resides in the regions of positively curved 
potential, d^V > 0. This occurs only in the region between two potential inflections around 
the minimum, which is broad enough to model the current epoch of cosmic acceleration if 
fa ~ 0{l)Mpi. However, for models where f a < 0(l)Mpi, there is no region of the convex 
part of the potential where the slow roll condition d^V <C V/Mpi holds. Even if the field 
starts at the inflection, the restoring force ~ d^V will be too large and pull it down very 
quickly, in time At inflation <, Hq 1 . Thus to mimic dark energy for longer such a pNGB 
had better start near a maximum. But in this region the pNGB is tachyonic. It grows 
exponentially, because its vev is sliding down the potential ridge. This tachyonic instability 
is slow if f a >j 0(l)Mpi, extending the longevity of the dark energy epoch. 

When f a < Oil) M pi the instability develops much faster. Generically the pNGB will 
not remain perched exactly on the top of the potential everywhere inside the horizon even 
if it had started there, because it will get displaced by quantum fluctuations during early 
cosmological inflation. The amplitude of these fluctuations sets the scale of the minimum 
initial displacement from the quintessence maximum, which will be of the order of the early 
inflation Hubble scale, < 10 13 GeV. For a given f a < 0(l)Mpi the field will roll down 
the slope fast. Suppressing this for long enough is what places a strong bound on f a of 
quintessence pNGB in this case. Further, the perturbations imprinted in the pNGB dark 
energy during early inflation will grow due to the tachyonic instability, and may lead to 
potentially observable large wavelength perturbations in the dark sector. In addition, since 
the pNGB vev will change in the course of cosmic evolution when f a Oil )Mpi, the dark 
energy equation of state will vary [16], possibly even leading to an increased rate of cosmic 
acceleration at late times [26, 27], and yielding increased supernovae dimming at small 
redshifts. This offers a novel possibility of correlating a variation of w(z) with the dark 
energy perturbations. 

Finally, we will also point out that pNGB dynamics could still accommodate the present 
epoch of cosmic acceleration even if f a < M P i - if there are several pNGB fields. The 
mechanism which lends to this is a scaled-down assisted inflation [28], whose high energy 
variant may emerge from string theory [29] . This can ameliorate some of the concerns about 
the possibility of embedding pNGB quintessence with f a ^> Mpi in fundamental theory 
[20, 21]. Instead of a single field pNGB quintessence, one may be able to use several fields 
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with f a < 0(l)Mpi to dominate the universe today. Such models may actually lead to 
a number of observational signatures that could be explored by the future probes of dark 
energy. 

Let us first briefly review the cosmology of pNGBs. The low energy dynamics of a pNGB 
is given by the Lagrangian (with the (— , +, +, +) metric convention) 

£ = -i(V0) 2 -/i 4 (l-cos(|))+.... (1) 

The ellipsis in Eq. (1) stands for the higher order corrections from the instantons as well as 
possible higher derivative terms. For now, we will ignore them. In an FRW geometry, the 
homogeneous pNGB mode obeys the field equation 

<j> + 3#0 + d^V = , (2) 

where H = a/a is the Hubble parameter and a the FRW scale factor. If the pNGB is in the 
slow roll regime, such that its potential energy dominates its kinetic energy, V > |(<90) 2 , and 
the restoring force in (2) is overwhelmed by the cosmological friction term, d^V < 9i? 2 /4, 
the cosmological dynamics is dominated by the potential, yielding cosmic acceleration. 

Naively this seems to occur when pNGB is displaced away from the minimum of the 
potential to ~ f a (modulo 2irf a ), such that H 2 ~ 0(1)/jl 4 /M pi (we use M Pi = 1/^8ttG n 
throughout) and d^V ~ ji A / ' f 2 , suggesting that as long as f a > 0(1) Mpi the slow roll regime 
may set in [13] . The scale of inflation driven by pNGB depends on the scale /i generated by 
non-perturbative symmetry breaking [13, 14, 19, 20, 23], and so to generate a sufficiently long 
epoch of cosmic acceleration today one must require ji ~ 10 _3 eV, such that tjiq < 10~ 33 eV, 
as mentioned above. 

However this argument is a bit too quick. Slow roll regime may only occur if the potential 
is convex [17, 18]. Hence the field should be within the interval around the minimum bounded 
by the two adjacent inflexion points. In this regime, the slow roll condition Mpi d^V/V <C 1 
gives 

Mpi sin (<j)/f a ) 

f a 1 - COS (<f>/f a ) ■ 1 ] 

If we restrict ourselves to the region for which V ((f)) is convex, the left hand side of the 
above equation attains its minimum at the inflexion points (^inflection — tt f a /2 ± nir f a . 
There, its value is Mpi/f a , so that one can hope to realize slow roll only if f a is larger than 
Mpi. Otherwise, the restoring force ~ d^V will pull the pNGB down to the minimum 
of the potential too quickly. Indeed: although the curvature of the potential in (1) is 
dlV = ^2 cos(0// a ), and it is dominated by friction near the inflections, the potential is well 

™ J a 

approximated by the linear term V ~ /i 4 (l — A(p/f a ) there. Soon after the field is released 
from an initial position near the inflection, its velocity saturates at <j) ~ —fj, 2 Mpi/(^/3f a ), 
so that the kinetic energy of the field is comparable to its potential unless f a > Mpi. This 
will happen after the vev of the field moves by ~ f a . Approximating the field acceleration 
by (j) ~ — fi 4 /f a , we can estimate the duration of an inflating phase in this case to be 
Ati nfla tion ~ jp, or therefore 

Atinflation ~ -TJ~ H Q 1 . (4) 

Mpi 
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Thus this regime will not last long enough unless /„ > M Pl . 

What happens if the pNGB field starts near the potential maxima? Clearly for as long as 
it sits precisely at the maximum, it would behave exactly like a cosmological constant. If such 
a state were very long lived, it may be able to imitate dark energy even if f a < M Pi . Getting 
it up there clearly exacerbates fine tuning. However one should first develop the means for 
quantifying the degree of such fine tunings. One might even be able to justify a special value 
of the pNGB vev dynamically: e.g. in the early universe before inflation the pNGB was a flat 
direction, and so it could have fluctuated a great deal. Its vev could have assumed any values, 
mostly away from the value required to mimic cosmological constant, = (ir + 2nn)f a for 
some integer n. However, the required field vevs may have been attained in some regions of 
the meta-universe simply because of continuity, and then frozen-in by early inflation. When 
the pNGB potential is generated later on, this region acquires a small vacuum energy in the 
pNGB sector. So we may just happen to live in one such domain - inside a pNGB domain 
wall - expanded to tremendous sizes by inflation. Similar mechanisms are also used as a 
catalyst of chaotic inflation [17] or topological inflation [31]. 

The problems with such arguments arise when one considers dynamics. If the initial vev 
starts inside an interval of width nf a centered at maxima at (n + 2nix)f a and bounded by 
the adjacent potential inflexions, the pNGB is a tachyon: near its maxima the curvature of 
the potential is negative. The pNGB wants to slide off of the peak of the potential, and its 
vev depends exponentially on time. When f a < M Ph this time variation is fast: the field 
changes exponentially in a Hubble time. Hence in this regime, the field is never in slow roll 
in the strict sense, and the epoch of dark energy domination will quickly cease. This is really 
the main reason that pNGB quintessence models work better with large f a , as we will now 
elaborate. 

To illustrate this, we can solve the field equation (2) by approximating the potential as 

v = { 2/i 4 - §m 2 (0 - (j) max ) 2 + . . . , for <p max - nfa/2 < < (j) max + nf a /2 ; 
~" \ \m 2 ((t) - (f) min ) 2 + ... , (f) min - TTfa/2 <(f)< (f) min + Txf a /2 . 

Here 4> ma x = (vr + 2nn)f a for some integer n, whereas (f> m in is a minimum adjacent to it, 
where pNGB wants to settle in, and where we have defined 

m2 = \d}y((f) extrema ) \ — — . (6) 

J a 

The corrections to the potential (5), denoted by the ellipses, are small away from the in- 
flections, but become important near them since they are smoothing the potential towards 
the exact cosine form. Nevertheless the corrections do not significantly alter the qualitative 
aspects of the conclusions and so we can ignore them for the most part. Now, if the vev 
is initially displaced from a maximum <f) m ax to 0j = (tt/2 + n7r)f a + <po, where ^ <po < fa 
so that dfrV < 0, and if the pNGB sector dominates the universe, the dynamics will be 

well approximated by + 3H<p — m 2 <j) = 0, with H = ^j2/3fi 2 /Mpi ~ const, with solutions 

oc exp(a±t) where ct± = [-3H ± \f^H 2 + 4m 2 ]/2. The leading order late time behavior 
is dominated by the growing, tachyonic mode oc exp(a + t). The power a + depends on the 
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ratio of H/m ~ f a /M P i, being 2 

for f a <0(l)M Pl ; 

f a >0(l)M Pl . [<) 

In either case, if the field had started at 0« with vanishing velocity, then at a time t >, l/a+ 
its vev will have moved to approximately 

<f>{t) ^<p m a X + \<Poe a+t + ... . (8) 

Once the total displacement at a time t becomes of order /„, A0(t) ~ \4>oe a+t ~ /„, the 
field will have approached the inflections of the potential, where the corrections omitted in 
(5) become important. 

The effective mass m^j = <9|V changes along the approach to the inflection, eventually 
vanishing there. Nevertheless the kinetic energy the field accumulated on the approach and 
the restoring force, oc d^V ^ there, will push the pNGB through the inflection, ensuring 
that it does not stop there. By this time the field velocity reaches ~ a + A(p(t) ~ a+/ a) 
and so the kinetic energy becomes <fi 2 ~ c^+fa- Depending on the ratio H/m ~ f a /Mpi the 
field may stop to behave as vacuum energy by this moment. Indeed, the equation of state 
parameter w(t), defined as w(t) = ^J^v during this epoch grows in time approximately as 

»W--l + i(^) 2 + -, 0) 

so that by the time A0(t) ~ f a , the parameter w becomes 

Wend ^ -I + J — ) +••• • 10 

Note that if m/a + ~ 0(1), or equivalently f a < G(l)Mpi (7), w will move significantly up 
from —1 by this time, and in this case cosmic acceleration will cease as the pNGB arrives at 
the inflection. 

After the inflection, the field dynamics is well described by <p + 3H<fi + m 2 = 0. What 
happens to it depends sensitively on the magnitude f a . If /„ >, Mpi so that m < H, 
from (10) it is clear that by the time the field passes through the inflection the universe 
is still accelerating (albeit marginally for m ~ H). By virtue of f a >, M Pi the field next 
encounters a sufficiently broad swath of convex and almost flat potential plateau, where 
it will in fact dissipate at least some of its kinetic energy acquired during the initial fast 
roll stage, somewhat prolonging the phase of acceleration. In the converse case f a < Mpi, 
Eq. (10) shows that the cosmic acceleration has all but stopped. In this case the potential 
plateau which follows the inflection is too narrow to slow the field down. Due to its kinetic 
energy, comparable to the potential at inflection, the field will zoom by, yielding only a 
small fraction of an e-fold. In either case, eventually the field will begin to oscillate around 
the minimum of (5), obeying a well known asymptotic attractor solution given by H ~ 



_2_ 

3t' 

2 Notice that the slow roll approximation would yield <j> oc exp(at), with a = m 2 /3iJ, that corresponds 
to a good approximation of the correct result (7) only in the case f a > 0(l)Mpi. 
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4> — 4>min + Acos(mt + 5) /a 3 / 2 . By virtue of the virial theorem, the field energy density in 
this epoch scales as p p ngb — /-^(^fy) 3 ) where a is roughly the value of the scale factor when 
the epoch of acceleration ended. Depending on the scales however, it may yet occur that 
after the first pass through the minimum the field continues to climb up the potential past 
the inflection on the opposite side. In that case, the field would yield another short burst of 
inflation as it slowly climbs the opposite slope, before returning to oscillations around the 
minimum. The possibility of having inflation briefly start again after the passage through 
the minimum has been considered earlier in relation to the cosmological moduli problem 
[32]. Such dynamics could simulate the behavior of dark energy with w < — 1, as recently 
discussed in [26, 27]. 

To summarize this discussion, we see that a universe controlled by a pNGB may first 
undergo a stage of inflation during which the field changes quickly, akin to the fast-roll 
inflation of [33], which may (for f a > Mp{) or may not (for f a < 0(l)Mpi) be followed by 
a stage of slow roll inflation, and then eventually end up like a universe dominated by an 
oscillating massive scalar field. If m/a + ~ 0(1), or therefore f a < 0(l)M P i (7), cosmic 
acceleration will come to a close by the time the pNGB vev gets to the inflection. In this 
case the duration of the accelerating epoch is fully controlled by the coefficient a+ ~ m (see 
Eq. (7)) and the initial displacement from the maximum, 

A*~— ln(^). (11) 

m v 0o ' 

If a single pNGB is to mimic dark energy today, At must be of the order of the current age 
of the universe, ~ l/H . To ensure this the initial displacement of the pNGB vev away from 
the potential maximum should be no more than 

0o max ~ f a e- m/Ho ~ f a e~ MM , (12) 

to ensure that the universe is accelerating now. 

Herein lies the problem: the condition 0o < 0o max will not hold unless f a >, Mp\. By 
its very definition, pNGB couplings are very weak and its potential, generated by non- 
perturbative effects, is very flat 3 , compared to a typical scale of early inflation. Therefore 
during early inflation which shaped our universe and generated primordial perturbations, the 
pNGB field is relativistic, with d^V <C Hj. Hence like any relativistic field during inflation 
it must have fluctuated a lot around its mean value. Its dynamics was stochastic [18, 34]: 
imagine that in some region of the early universe the pNGB vev did attain exactly the value 
corresponding to the maximum of the potential in (1), = (n + 2nn)f a , when inflation 
started. Because it is dynamical and relativistic, the field subsequently gets displaced by the 
quantum fluctuations generated in de Sitter environment [18, 34, 35, 36]. It "wiggled" a little 
inside the inflating region from one Hubble time to another. The scale of these quantum 
fluctuations is given by the root mean square value of the field around its expectation value, 
set by the two point function of a relativistic quantum field in de Sitter space [18, 34, 36]: 

<0 2 t >= (-1) / — . 13 

v 2n ' JH I e- H i t k 

3 In fact at energies above the scale of shift symmetry breaking, there is no potential for pNGB at all. 
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From this, the jump of the field between two subsequent Hubble intervals is 



ni 

r 2tt 



(14) 



The constraints on f a when it is smaller than 0(l)Mpi, derive from this. While dependent 
on the scale of early inflation, they are nevertheless useful. If the field is initially displaced 
by ~ Hi/2% (14), due to the tachyonic instability it will slide off the maximum in time At 
given by (11). If this is to yield enough late cosmic acceleration, At >, 1/H . Thus we find 



where we used (6) and H ~ fi 2 /Mpi. The inequality (15) implies that the logarithm ranges 
between ~ 15 and ~ 36 for the range of Hj from about 10 13 GeV down to about 10 4 GeV ', 
which cover many generic inflationary models including chaotic inflation [17] and hybrid 
inflation [37], its supergravity realizations [38] and modern hybrid inflation models on the 
string landscape [39]. Using (15) we can estimate that f a must be greater than about a tenth 
of Mpi if pNGB is to mimic dark energy today: 



In fact, numerical considerations show that the bound is even stronger by a factor of a few. 

In Figure 1 we plot the current value of the equation of state parameter of the pNGB, 
as a function of its initial value O an d of f a /M P i. Current observations relegate the allowed 
part of parameter space to the top right corner of the figure, and - barring fine tuned initial 
conditions - a typical bound on f a reads f a > 0.2 Mpi, to fit the observations. We should 
remark, here, that the boundary between the forbidden light-colored region and the allowed 
darker one is not as sharp as it appears on the plot. Indeed, for values of 0o and f a close 
to the boundary, we have regions where the pNGB has performed a half oscillation around 
the minimum before the present time, and has been recently climbing its potential 4 . In such 
cases, the initial potential energy of the field starting near the maximum does not dissipate as 
quickly, and they swing past the minimum beyond the opposite inflection, where they linger 
a little on the potential. The study of the viability of these regions is more complicated, as 
the full time dependence of the parameter w(a) should be taken into account in this case, 
along the lines of [40]. These regions could actually provide a better fit to supernova data, 
mimicking w < — 1 [26, 27]. 

We can thus conclude that models where f a < 0.1 Mp\ are extremely fine-tuned: to 
get a pNGB to mimic dark energy when f a fails to satisfy (15) one needs extreme fine 
tunings of the initial state of the pNGB. The field must start in a special state where the 
subsequent quantum fluctuations precisely cancel against the prearranged initial dark energy 
inhomogeneities. This is extremely unlikely, and goes against the grain of the standard 

4 These regions are not represented in Figure 1, because in the plot we have computed the value of w at 
the moment in which f2 ma tter crosses the value 0.3 for the first time. In the "forbidden" region, the value 
of w is larger than —0.7 when this first crossing happens. But for values close to the boundary between the 
allowed and the forbidden region, Q ma tter actually performs small oscillations around the value 0.3, while 
the equation of state parameter w of the pNGB oscillates between —1 and 1. 




(15) 



(16) 
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inflationary lore. In inflation, such special features of the initial state would be erased 
by exponential expansion, and the subsequent structure would be imprinted by quantum 
fluctuations, leading to (15) for pNGBs with f a < Mpi. 

If we parameterize the end of inflation by picking a spacelike surface along which the 
adiabatic perturbations generated during inflation are 5p/p ~ 10~ 5 , we will find that the 
distribution of the vevs of the pNGB field is also slightly inhomogeneous, given by an approx- 
imately scale invariant spectrum at scales >, 1/ Hi and with the field perturbation estimated 
by, roughly, Eq. (14). These perturbations of the pNGB vev will remain frozen until the 
horizon reentry at a much later time. They will yield small initial perturbations of the dark 
energy density, 

s _pm^w~J?L (17) 

Pde V 8^' 1 ) 

where we assume that the field was initially exactly on top of the potential. In evaluating 
the last line of (17) we have used (6). If <fi started off farther from the maximum, the 
perturbations would be larger, 

Pde ^Ja 

Initially, the dark energy perturbation may be very small. For a typical inflationary model 
with Hi < 10 13 GeV, if the pNGB decay width is large, f a >, Mpi, the initial dark energy 
perturbation would be Spde/pde ~ 10~ 14 — 10~ 7 , subsequently slowly growing in time, as 
the pNGB slides off the maximum. When f a < 0{l)M Ph however, the perturbation would 
both be initially larger and would grow faster. Indeed, in regions of the postinflationary 
universe where inflation displaced pNGB away from the maximum by ~ Hi the late stage 
of acceleration driven by the pNGB would terminate faster than in those regions where the 
field remained perched on top of the potential. Such models can develop large perturbations 
at late times in the dark sector, because the initially small dark energy density perturbations 
will grow as pNGB slides down the potential 5 . Even if they start as small as (17) or (18), 
by the time the field approaches the potential inflection, A<fi — > f a , the perturbations will 
increase according to 

6PDE e "^( A0(At) ) 2 ^0(l), (19) 

PDE end V fa ' 

increasing as f a drops below Mp\. Thus the dark energy perturbations could in principle 
grow to become sufficiently large today to compete with the inflationary perturbations at 
large scales. We expect that these dark energy perturbations would yield isocurvature per- 
turbations at very large scales and very late times, and the mechanism would be akin to the 
curvaton mechanism of [42]. While here we will refrain from a detailed computation of the 
spectrum and evolution of dark energy perturbations, we note that such a study has been 
performed on pre-WMAP data in [43] (a previous study of this issue can be found in [44], 
where the authors, however, focused only on the special value f a = M P {). The constraints 
found in [43] are weaker than the ones imposed by supernovae. The reason, roughly, is that 
during inflation structure grows slowly, and hence the perturbations remain frozen until the 
end of the accelerating epoch. An analysis that includes the more accurate WMAP data 



5 For some choice of parameters, the fluctuations might backreact on the evolution of the zero mode of 
the pNGB [41], affecting its equation of state. 
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may help reducing the available portion of the parameter space, 6 although it should also be 
clear that if f a >, 0(1) Mpi such effects may remain small 7 . The possibility of generating 
and observing the imprints of early inflation in dark energy has been recently discussed in 
more general terms in [47]. 

These arguments show that the pNGB quintessence models need f a > 0(l)M Pl to be 
viable. How could such models arise from fundamental theory, given the problems so far 
with deriving large f a in string theory [20, 21]? In fact, there is a simple way to provide 
for this using pNGBs with f a < 0(l)Mpi, if there is many of them 8 . It is based on the 
mechanism of assisted inflation [28, 29]. Imagine that there are M copies of a pNGB with 
f a < 0(l)Mpi, which interact with each other via gravity. Then the total Lagrangian of the 
system is the sum of M copies of (1). Assume for simplicity that all the fields are displaced 
away from the minimum by the same amount. In that case, we can introduce a collective 
coordinate 

^ = ^0, (20) 

whose effective potential is given by 



K//=AV(l-cos(-^r)). (21) 



Thus we see that the collective coordinate tp behaves as a single pNGB with a potential 
controlled by the symmetry breaking scale /i e // and decay width f e ff 

toff = -A^V, 

feff = AT 1 ' 2 fa- (22) 

Therefore the total vacuum energy and the effective mass analogous to (6) for the collective 
mode tp are 

X ef f = J\ffl 4 =J\f\, 

2 A/> 4 _ m 2 

m eff ~ Jpj2 ~ ~Jf ■ W 

We can therefore approximate the potential for tp with an expression just like (5), where we 
replace fi, m and f a by /i e //, m eff and f e ff given in (22), (23). The analysis of the cosmo- 
logical dynamics controlled by the collective mode tp then is the same as the investigation 
of the cosmology dominated by a zero mode of a single field <f>, and we see that by virtue of 
Eqs. (22), (23) the duration of the epochs of dark energy is prolonged. If the field tp starts 
below the inflection, where the potential is convex, the duration of the slow roll regime is 
approximately the same one that we would find in a simple quadratic potential. Indeed, the 
field tp will be slowly rolling until its vev is at a distance of the order of the Planck mass 

6 See [45, 46] for an interesting discussion of the possibility that such perturbations might play a role in 
explaining the observed lack of power in low-£ CMB multipoles. 

7 When f a £ 0(l)Mpi the pNGB vev will change more slowly in a Hubble time ~ 1/Ho, and also the 
field need not start from the maximum of the potential. 

8 Wc exploit the possibility of having many pNGBs with couplings which differ from each other, in contrast 
to the case of two pNGBs of [22] where an accidental symmetry relating axion couplings was used to generate 
a large f eff . 
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from the minimum of V (</?). The total number of e-folds will be ~ flff/M Pl ~ A/", which 
we can again estimate using the formula relating the number of e-folds N and the field vev 
for a polynomial potential V ~ \(p n /n, N = (y? 2 — ip^) / (2nM Pl ) [17, 18, 30]. Assuming 
J\f > (Mpi/ fa) 2 we can approximate the potential by the quadratic term, as opposed to the 
linear potential as in the case leading to Eq. (4). This yields 

Notice that this is the regime considered, in the context of primordial inflation, in [29]. 
Likewise, the tachyonic instability is also slowed down. When Af > MpJ fl so that f e ff > 
Mpi, the value of a + in Eq. (7) is now given by a + ~ m^ff/SH, and so the Eq. (11) is now 
replaced by 

A ^ifMn(^UA^ln(^). (25) 
mi ff V V?o / m 

Again, it takes f 2 ff/M Pl e-folds to get to the inflection. But now even if f a /Mpi < 1, as long 
as there are many pNGBs such that VJ7f a /M P i > 1, a full e-fold of acceleration or more 
may ensue. 

Obviously it is important that the pNGBs are not directly mutually interacting. Direct 
interactions are detrimental to the effect of assisted inflation, which can be seen very simply 
from the following argument. The effect of assisting works by taking a given potential of 
fixed curvature and then prolonging the epoch of dark energy supported by this potential by 
taking J\f degenerate copies [28]. If the direct interactions are then turned on, they will lift 
the degeneracies of the mass spectrum of the fields. Some fields will be lighter, but most will 
get heavier. This is akin to a generic effect of mass splitting of a degenerate quantum system 
under a perturbation. Hence, in spite of a large number of fields, some of them will become 
so heavy under the interactions that they will not be able to stay on the potential ridge, 
and will tumble down. However, it may be possible to ensure that the mutual interactions 
are miniscule by arranging for the right compactification of string theory, where pNGBs are 
axion-like degrees of freedom present in the moduli space, along the lines of [29]. A version 
of such a mechanism with very low mass scales may find its home in the landscape of string 
solutions. 

Our argument that many copies of pNGBs prolong the epoch of dark energy relative to 
a single one is clearly sensitive on the initial conditions for the individual pNGBs in the 
early universe. The assisted dynamics works best if the pNGBs start out at roughly the 
same place, away from the minimum of the potential, where they hold each other up on 
the potential slope. In general, however, the effective field cp will have some width, since 
the pNGBs may start at a different place. Further, they may also have slightly different 
effective masses (even if their \x and f a are all the same, the effective mass is determined by 
the curvature of the potential and hence depends on the position of the field). Because of 
this, the field f will eventually begin to fall apart in the course of evolution, and this will 
shorten the dark energy epoch and may lead to a new source of dark energy perturbations. 
While it is obvious that tuning the initial values of the separate pNGBs will hold this under 
control, it would be interesting to revisit this issue in more detail. 
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Finally, let us comment on the idea that the pNGBs might have so extended a spectrum 
of masses and couplings that some of them contribute to dark energy today while others have 
already been oscillating for a long time around the minimum of their potential, behaving in 
effect as dark matter. In this scenario, dark matter and dark energy would share a common 
origin, a fact that might help with the coincidence problem. We find that such mechanisms 
are hard to implement naturally in a simple setting. To illustrate these difficulties, consider 
the simplest case of two pNGBs and ask that one of them (with parameters fi\ and j\) be 
dark matter, while the other (//2, $2) be dark energy. The first pNGB will start rolling along 
its potential when the Hubble parameter ~ T 2 /Mpi equals its mass fif/ fi- This happens at 



a temperature T\ ~ ji\J Mpij f\. The pNGB subsequently evolves as matter, so that 



and it starts dominating the energy density of the Universe at the temperature T e ~ 1 eV 
of matter-raidation equality. We thus find the relation /ii ~ T e {Mp\j fi) 3 ^ 2 . By requiring 
T\ > T e , we get f\ < Mpi, that implies fii > T e . The parameters of the second pNGB are 
determined by the requirement the it is dark energy today: ji2 — 10~ 4 T e , f 2 — Mp\. To 
solve the coincidence problem, one would need to relate the two pNGBs in some natural 
way, requiring that their mass scales and couplings fik, fk are related by a common new 
physics at some new scale A. Parameterizing this relationship by a simple power law, so 
that /ii j2 = A 1_/3 /^ 2 , we get the relation 



As fi < Mpi, (3 has to be negative, and [i should be thus inversely proportional to some 
power of /. Typically models of pNGB have the opposite behavior, and so this would require 
some complicated and unusual construction to reproduce the desired spectrum of pNGBs, 
pursuit of which is beyond the present work. 

To recapitulate: we have revisited models where dark energy is dynamical, controlled by 
a pNGB field. We have stressed that unless f a >, 0(l)Mpi, models with a single, or very few, 
pNGB(s) are not good explanations of the present epoch of dark energy domination because 
of the tachyonic instability plaguing most of the field initial values. Even if the field starts out 
very close to a maximum of its potential, because it is so light quantum fluctuations during 
early inflation will shift it away from the maximum by an amount proportional to the Hubble 
scale of inflation Hj. This would exacerbate the fine tuning problems of dark energy. Given 
the concerns raised about embedding pNGB models with f a > M Pi in fundamental theory 
[20, 21] on the one hand, and the acute absence of technically natural models of quintessence 
on the other [20, 12], this may present a challenge for building dark energy models within 
the realm of effective field theory consistent with fundamental theory. We have noted that 
a possible solution of this problem can be developed by scaling down the mechanism of 
assisted inflation [28, 29] to very low scales. If there are several ultralight pNGBs, they may 
drive very late epoch of cosmic acceleration together even if their individual decay constants 
satisfy f a < M P [. We should note that although the assisted inflation idea does not seem 
to help other quintessence models [48], it works well for pNGBs. If quintessence is really 
an admixture of many modes, its decay may leave signatures in the form of inhomogeneities 





(26) 



{h/M Pl f 2+ P = 10 4 . 



(27) 
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of dark energy distributions. These inhomogeneities are triggered by quantum fluctuations 
during early inflation, just as the primordial density perturbations, and ultimately baby 
universes in self- reproducing inflationary models [34] . Thus in some sense such perturbations 
may carry information about the cosmic seeds, frozen in the cosmic permafrost that is dark 
energy. It would be interesting to carry out more precise calculations of the amplitude 
of these perturbations and check their observational consequences, since they may affect 
cosmological anisotropies at the largest of scales [45, 46]. Further, in such models dark 
energy dynamics may also result in a variable equation of state w(z), which could account for 
additional dimming of supernovae at low redshifts [26, 27], warranting coordinated searches 
for w' and for dark energy perturbations. 

Finally, we cannot refrain from reminiscing about one other question, often heard in the 
field of dark energy: namely, "Who ordered this?". The pNGB quintessence models strive 
to be natural in the sense of effective field theory, since the scales and couplings governing 
their dynamics are radiatively stable. However the pNGB dark energy with f a < 0(l)Mpi, 
which were a central subject of this study, do not appear all that natural or simple from 
the point of view Occam's razor, a criterion often employed to justify or belittle a model of 
dark energy. The oft-quoted versions of Occam's razor are commonly extrapolated to imply 
that dark energy may be the simplest of possibilities: a cosmological constant. Yet this is 
a heuristic argument that does not necessarily give correct answers. It is merely a loose 
guide enabling one to pick a scientific hypothesis based on the fewest currently unproven 
assumptions. To fully accept arguments based on simplicity and Occam's razor, we need 
to go beyond naive and develop a complete theory of the dark sector, that goes beyond 
the taxonomy of cosmic coincidences [10]. Before we have it, we may be dismissing ideas 
too quickly. The appearance of simplicity may be deceiving, as witnessed by a dearth of 
complete theories of vacuum energy. Just putting a number A in the equations is merely 
a fit. The pNGB quiA/tessence models are motivated by an attempt to simultaneously 
address two theoretical issues: (1) protect the ultra-low scales of quintessence from radiative 
corrections involving the rest of the world 9 , and (2) derive ultra-weakly coupled pNBG-like 
fields, simulating f e ff > Mpi, from fundamental theory. If it so happens that such models 
bring about a number of observational consequences worth studying, all the better. Similar 
opinions about Occam's razor have been expressed before on numerous occasions (see [49] 
for a nice summary). We end by quoting one such example: 

"There are more things in heaven and earth, Horatio, 
Than are dreamt of in your philosophy. " 

William Shakespeare, Hamlet 

Whether or not it applies to our world remains to be seen. 
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Figure 1: Current value of the equation of state parameter wq of the pNGB. Going inwards 
from the top right corner of the figure, the different colors correspond to areas with — 1 < 
w < -0.99 (red), -0.9 < w < -0.99 (orange), -0.7 < w < -0.9 (green) and w > -0.7 
(yellow). 
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